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Complete ^-matrix in a microwave cavity at room temperature 
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We experimentally study the widths of resonances in a two-dimensional microwave cavity at room 
temperature. By developing a model for the coupling antennas, we are able to discriminate their 
contribution from those of ohmic losses to the broadening of resonances. Concerning ohmic losses, 
we experimentally put to evidence two mechanisms: damping along propagation and absorption 
at the contour, the latter being responsible for variations of widths from mode to mode due to its 
dependence on the spatial distribution of the field at the contour. A theory, based on an S-matrix 
formalism, is given for these variations. It is successfully validated through measurements of several 
hundreds of resonances in a rectangular cavity. 

PACS numbers: 05.45.Mt, 05.60.Gg, 05.40.-a 



I. INTRODUCTION 

In the field of Quantum Chaos, microwave experiments 
have proved to yield very important breakthroughs in 
providing versatile analog models of quantum systems in 
the domain of classical electromagnetic waves yj. Room 
temperature experiments have opened the way[2| rapidly 
followed by experiments in superconducting cavities |3j. 
In a first stage, studies have mainly been concerned with 
the verification of predictions issued from Random Ma- 
trix Theory or from semiclassical approaches regarding 
spectral fluctuations. Losses, which were originally ab- 
sent from theoretical models, were seen as severe draw- 
backs in the seminal experiments, especially for an accu- 
rate analysis of resonance frequencies (see e.g. 0,0,0)- 
The first account of resonance widths observed in su- 
perconducting cavities was related to coupling losses in 
the absence of ohmic losses and measuring widths essen- 
tially amounted to measuring intensities at the locations 
of few antennas 5j. During the last decade, the great 
flexibility of microwave cavities has led to an important 
diversification of geometries and configurations in order 
to investigate the spectral correlations and the spatial 
distribution of the field, in closed or open, disordered 
and/or chaotic cavities (see pj for a review). Neverthe- 
less, until recent years, the impact of the different loss 
mechanisms, present in these systems, on their spatial 
or spectral statistical properties had obtained very lit- 
tle consideration. Indeed, as long as losses are weak, 
resonances can be viewed as isolated. On the contrary, 
for increasing damping, resonances are no longer easily 
distinguished due to modal overlapping and the very de- 
scription of the wave system in terms of modes loses its 
pertinence. Since the seminal papers by Ericson in nu- 
clear physics and by Schroeder in room acousticsQ, 
the regime of large modal overlap has been abundantly 
studied in the context of quantum chaos[3. l9l [Toj . 

The question of intermediate modal overlap for which 
resonances can be distinguished but broadening is no 
longer negligible is essentially open as yet (see the ex- 



cellent review |lljh In the present paper we propose to 
help pave the way of a more complete understanding of 
microwave cavities at room temperature by accounting 
for the presence of essentially two kinds of loss mech- 
anisms, namely ohmic damping at the boundaries and 
coupling to the outside through antennas. To be able to 
separate their respective contributions to the broadening 
of resonances, a thorough analysis is required of the way 
the wavefunctions are spatially distributed throughout 
the cavity. 

The cavity we have actually used for our experiments 
is composed of two rectangular OFHC copper plates be- 
tween which a copper rectangular frame is sandwiched. 
The rectangular frame has been machined as one piece 
and serves as the contour of the cavity. The cavity 
may thus be viewed as the slice of a rectangular waveg- 
uide closed at both ends, with contour C of length L = 
2.446 m, section S of area A = 0.3528 m 2 and thickness 
d = 5 mm. As long as the wavelength A is larger than 
d, the boundary conditions in the z direction (perpen- 
dicular to the top and bottom plates) only admit Trans- 
verse Magnetic (TM) two-dimensional (2D) modes. The 
whole structure is tightly screwed and 10 holes have been 
drilled through one of the plates to introduce 10 anten- 
nas, which protrude a length I into the cavity. The an- 
tennas are monopolar with SMA connectors which are 
commonly used in the frequency range from to 18 GHz. 
The positions of the antennas are displayed on figure Q 
For a measurement, only one antenna at a time is used 
as a microwave emitter and another (in transmission) or 
the same (in reflection) as a receiver. The other unused 
antennas are terminated by 50 O loads so that all anten- 
nas behave the same way regarding the losses they imply. 
These antennas are linked to an HP 8720 D vector ana- 
lyzer through flexible cables. All the measurements are 
performed after a proper calibration to get rid of any par- 
asitic influence of cables and connectors and even of the 
analyzer itself. The measurements are given in terms of 
scattering coefficients which form the S-matrix ( g** s)> 
where Six (resp. S22) measures the reflection on port 1 
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FIG. 1: Schematic view of the microwave cavity with locations 
of antennas. 



the time evolution being conventionally written 
exp(— iu)t). 



Cavity 




-l* + l 



(resp. 2) and S12 (resp. S21) measures the transmission 
from port 2 (resp. 1) to port 1 (resp. 2). 

In the following section we develop an S'-matrix for- 
mulation for an ideal cavity through the introduction of 
an electromagnetic model of antennas which enables us 
to write the response of the cavity within the form of a 
Breit-Wigner decomposition. Then, in the third section, 
we complete this description with a perturbativc evalua- 
tion of ohmic losses at the walls of the cavity. We show 
that the resulting ohmic width of each resonance may be 
decomposed in two qualitatively and quantitavely dis- 
tinct contributions, one of them being sensitive to the 
spatial distribution of the wavefunction at the contour. 
Then, in section HVI we proceed to an experimental val- 
idation of our model in the case of a rectangular cavity. 
We show that, for each measured resonance, we are able 
to discriminate quantitatively among the two ohmic con- 
tributions to the total widths and the contribution due 
to the presence of the antennas. 



II. S-MATRIX FORMULATION FOR A CAVITY 
WITHOUT OHMIC LOSSES 

A. Electromagnetic model of antennas 

As an antenna, we use the terminal part of the center 
conductor of a coax (see Figure |2J). Far from this termi- 
nation, in a coaxial line, only transverse electro-magnetic 
(TEM) modes can propagate and the field results as the 
superposition of incoming and outgoing parts. In the 
vicinity of the termination of the line, hereafter called 
the perturbed region, perturbativc non-propagating waves 
exist 01 ■ The longitudinal variable z along the line is ori- 
ented outward from the cavity and its origin located at 
the border between the TEM and the perturbed regions, 
i.e. at a distance I* from the end of the antenna (see Fig- 
ure^. In the perturbed region, assuming a sinusoidal 
behavior, we write the stationary current I per t(z) as 

WO) = 1+^" + l-c~ ikz for - I* < z < , (1) 



FIG. 2: Schematic cut view of a coupling antenna. Physical 
regions introduced in our model are displayed along with the 
associated characteristic lengths. 



The current I(z) in the TEM region reads : 



Hz) = ■£ {Aout^ z +A in e- 1 n for z > 0, 



(2) 



where Z is the characteristic impedance of the coaxial 
line. The following continuity conditions are then im- 
posed 



Ipert(—l) — 

W(o)=/(o) 



Therefore 



I pert (z) = 1(0) 



sinfc(r + z) 
sinfcZ* 



pour 



r < z < 



with: 



7(0) = -|(A»t+ Ai n ) 



(3) 



(4) 



(5) 



Inside the cavity, time-independent Maxwell's equa- 
tions yield the following wave equation for the electric 
field E: 



AE + k E = iui/iJ . 



(6) 



where, assuming a pointlikc antenna at location f*o, the 
current density within the plane of the cavity reads: 

J(z,r) = I{z)5(r-r )z for -l*<z<0. (7) 

The 2D formulation of our problem is obtained by in- 
tegrating equation © along z in different ways for the 
left-hand side and the right-hand side. Indeed, while the 
left-hand side is easily integrated over the thickness d 
of the cavity, i.e. for — I* + I — d < z < —I* + I, the 
integration of the right-hand side is more involved. To 
account for the effective coupling of the electric field with 
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the current in the perturbed region, we define a coupling 
function f(z) on the interval [—1*, 0] which multiplies the 
current before integrating. To our knowledge, only nu- 
merical approaches of this problem have been published: 
using FDTD0] or modal decomposition 0, 0|. Here 



we adopt an effective description by assuming f(z) = 1 
on the interval [—1*, —I* + l e ff] an d fi z ) = on the rest 
of the interval. The length l e ff is an adjustable parame- 
ter lying between I and I*, most likely close to I. Hence, 
for M identical antennas, equation JfjJ becomes 



M 



iV Z sin z 



(A + k 2 )E z {r) x d = iu,/x £ 5(r- r c )^(A c out + A c in ) I 

2 fcief/ M 



-r+i. 



ft sinfc(P +z) 



sin kl* 



Z sin ^ cos ^ 

2 2 c— 1 



(8) 



where Zo = y/io/eo is the vacuum impedance. 

B. Breit— Wigner decomposition 

We now follow a standard approach in scattering the- 
ory to analytically express the ^-matrix of a microwave 
cavity coupled to pointlike antennas. This kind of cal- 
culations was initiated in nuclear physics |l7j and has 
been reproduced in various contexts since then (see e.g. 
references [IH, [H and [H|). The cavity is described 
as a closed system coupled to M channels, one for each 
antenna. The complete Hilbert space of the system com- 
prises the cavity and the channels. It is therefore decom- 
posed as the direct sum of Hilbert spaces associated to 
the inside and the outside of the cavity: £ = £i n © £ ou t- 



Though the inside Hilbert space is of infinite dimension, 
we adopt the commonly used simplification of a finite 
dimension N ^> 1 jnj. The Hamiltonian of the cavity 
TCi n is thus represented by an N x N matrix H . The 
cigenstatcs associated to Tti n are denoted The out- 
side Hilbert space is associated to the M antennas and is 
written as the direct sum: £ out {E) = £ i • -(B£m{E), 
where E is the energy of a continuum of scattering states 
denoted by \c,E) for channel c. Finally W denotes the 
coupling matrix of dimension N x M between the bound 
states of the cavity and the scattering states of the anten- 
nas. As long as the wavelength remains smaller than the 
distances between antennas, the direct coupling between 
channels may be neglected. The complete Hamiltonian 
H thus reads: 



N M M N , s 

H= / dE\c,E)E(c,E\ + J2J2(\fj,) / &EW» c {E){c, E\ + h.c. ) 



(9) 



w, 



out^-in ~t~ *out 



r 



The normalization conditions are: 

{u\n) = 5 V „ and (a, E\b, E') = S ab 5{E - E') . (10) 

The space representation of TLi n is given by: 

(r\H m \r') = -5(r-r')A P . (11) 

Likewise, the space representation of TL ut in the coax 
reads: 



(z a \Hout\Zb) = -(z a \Zb)- 



(12) 



For pointlike antennas the coupling will be represented 
by: 



(r\W out -, m (E)\z c ) = t c (E,z c ) 6(f- r c ) . (13) 



Let $ = ( va(*0 v/3i(zi) — v/3 M (zm)) denote a state 
of the complete system. Then the eigenvalue problem 
= Eq& may be written as: 
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A^ a (r) + J2c=i <H^ _ r c ) J dz c t c (E , z c )(pp c (z c ) = E tp a (r) 
t*[Eo, zi)tp a (fi) -(p'faizi) = Eoipp^zi) 

t* M (E ,z M )(p a (rM)-<Pp u (zM) = E (P/3 m (zm) 
I 



(14) 



where the prime symbols stands for the ordinary deriva- 
tive. 

According to the electromagnetic description given at 
the beginning of this section, t c (Eo, z c ) will be vanishing 
except in a perturbed region of length I* from the termi- 
nation of the antenna. It is therefore quite natural to fix 
the origin of z c at the limit of this region. For z c > 0, 
the field is written as the superposition of ingoing and 
outgoing waves 



1 



V/S.(*c) = -^( A o> 



■A c in G~ ikZc ) for z c > 0, (15) 



where k 2 = Eq. The factor 1/yk is required for dimen- 
sional and normalization purposes, and A\ n and A c out are 
dimensionless complex amplitudes. For z c £ [— Z*,0], the 
perturbed field 4>(z c ) in the coupling region is still to be 
a superposition of plane waves. Continuity conditions in 



agreement with our electromagnetic description imply: 

(16) 



0/J.(O) =V/J e (0) 



One therefore deduces 

/ » , » sin k(z r + 1*) . , . 

fo c (z c =WJ . , y- ' for - I* < z c < . 17 

sin ret* 

The M X M S'-matrix being defined by: 

Aout = StAin (18) 

where Ai n = (a]„ ... A™ ) T and A ou t = (-4Lt ■•• A IL Y ', 
only the functions ip/3 c (z c ) for z c > are relevant. Equa- 
tion HI 4(1 can thus be obviously reduced to: 



A(p a (r) +J2c=i S (* f ~rc) J^dZc t c (E ,z c )(f>p c (z c ) = E ip a (r) 



-V$ M (z M ) = E <pp u (z M ): 



(19) 



r 



where t c (Eo,z c ) = has been used for z c > 0. Here 
again we introduce the simplification of a non-vanishing 
constant value t c (Eo) for t c (Eo,z c ) only on the interval 
[— Z*,— I* + l e ff]- With this assumption and expression 
(ITTjl . the first equation in l(T§|) reads: 



M 



-A( Pa (^ + J2T c (k)5(r-r c )<pp a (0)=EoVa( i ?), (20) 



with 



c=l 



T c (k) 



t c (E ) sin 2 ^ 
k sin ^ cos ^ 



(21) 



Assuming identical antennas, i.e. T c {k) = T(fc), and by 
identifying JSJl with (J2DJ, one deduces 



T{k)=iVk^ _.-,± kl 



■ 2 kl e ft 
:!ir» w J 



Z sm ^ cos ■ 2 



(22) 



In spite of the reduction performed in equations 119(1 . 
by eliminating the source terms associated to the field 
inside the cavity flsLITfij . the condition of self-adjointness 
can be recovered through appropriate boundary condi- 
tions. In appendix IA II we show that it can be done 
through the following boundary condition for channel c: 



T*(fc)^(f c ) = ^ c (0) 



(23) 



We are now in a position to derive an explicit expres- 
sion for the S-matrix. If the energy dependence of the 
coupling is small on a scale of the order of the mean en- 
ergy spacing between neighboring modes, the S'-matrix 
can be written(Tll|: 

Sab = 5 ab -2iti{a,E\W\E~H e}s )- l W\b,E) , (24) 

where 

H eff = H -mWW 1 . (25) 
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The S-matrix can thus be rewritten: 

1 - iK 



S = 



1 + iK ' 



with 



K 



E -H 



:W. 



(26) 



(27) 



Standard linear algebra (see e.g. reference [Tg) trans- 
forms expression (gSJl into 



S = I M - 2mW 1 



1 



k 2 — H + mWW^ 



W . 



(28) 



Then, assuming a weak coupling, a pcrturbativc expan- 
sion to leading order (considering the isolated cavity as 
the "zeroth order" ) yields the following expression for an 
element of the S'-matrix: 



JV 



Sab = Sab - 2i7T ^ 



w*w ub 



(29) 



where the sum runs over the eigenstates of the isolated 
cavity with energies k 2 . In appendix I A 21 it is shown that 
the coupling matrix elements, in the space representa- 
tion, are given by: 



surface. If the conductor is perfect, E is perpendicular 
to the surface, H is parallel and there is no dissipated 
power — in the following, parallel or perpendicular will 
be understood with respect to the surface of the conduc- 
tor. This ideal situation will correspond to the zeroth 
order of our description of the field near the surface of 
the actual conductor. For a finite conductivity er, one 
can compute the first order corrections for the fields fol- 
lowing the standard approach described for instance in 
reference po| . 

To first order, the perpendicular electric field and the 
parallel magnetic field outside the conductor remain un- 
modified. Using appropriate boundary conditions to- 
gether with Maxwell equations, it may be shown that 
nonvanishing parallel components of both electric and 
magnetic fields exist inside the cond uctor. T hese fields 
decrease as exp(— z/5), where S = ^2/ fi c auj is the skin 
depth (/i c being the magnetic permeability of the con- 
ductor and a its effective conductivity) , and only depend 
on the zeroth order parallel component of the magnetic 
field H^p at the surface of the conductor. By continuity, 
one deduces the existence of a small parallel component 
of the electric field just outside the conductor: 



E 



^ = J^(l-i)(,A<>) 



(33) 



W u 



f{k)%{Tc) 
\Zk~TT 



(30) 



where ip^{r) = {r\[i) is the eigenfunction associated to 
k 2 . One finally obtains the following explicit expression 
for the S'-matrix elements: 

|2 



S a b 



2i«x 



JV 

E 



(31) 



{k 2 -kl + i\f(k)\ 2 Eo=x\Mrc)\ 



III. 



PERTURB ATIVE EVALUATION OF 
OHMIC LOSSES 



In this section we present the results deduced from a 
standard first-order perturbation approach whose valid- 
ity is restricted to nondegenerate modes (see Jackson's 
textbook H3). The power dP dissipated (ohmic losses) 
by a wave with frequency oj within the surface element 
da of a conductor is given by the flux of the real part of 
the Poynting vector through this surface. By adopting 
Jackson's convention exp(— itut) for the time dependence 
of the field, one has 



(32) 



where n is the unit normal vector directed toward the in- 
terior of the conductor and E and H arc the fields at the 



Using (|3^|) one finds: 



dP fj, c uS || g (o) M2 



(34) 



To obtain the total power dissipated through ohmic losses 
within a cavity, a mere integration of equation l|34|) over 
the walls is required. 

In an ideal 2D cavity, the electromagnetic fields does 
not vary along z: 



o 



£ (o) H and ^ (0) 







Ei°\x,y) 



(35) 



Denoting ip^ = Ei°\x,y) 1 the time-independent 
Maxwell equations are reduced to a 2D Hclmholtz equa- 
tion: 



(V? + e^ 2 )^ (0) =0 



(36) 



where the transverse gradient operator Vt is associated 
to the (x, y) coordinates, and e and \i are respectively the 
permittivity and the permeability inside the cavity. On 
the contour, tp^ obeys Dirichlet conditions: 



(0) 



onC. 



(37) 



This yields a complete analogy with the free propagation 
of a quantum particle in a 2D infinite well. This type 
of system is commonly called a quantum billiard. For a 
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given mode, the ohmic dissipated power reads: 



where 



P 



2aS 



d£ / dz||riAi? (0) || 2 o 



tfda\\nAHW\\l 



ends 



(38) 



One defines the integrals I\ and I2 so that equation H38|) 
is re-written: 



P 



h+h 



2aSfj, 

Using Faraday's law, one obtains: 

h=^etfda\^f 



where £ is defined by: 
1 



JC 



d£\d n *p 



I^^JJ da^ 



(0)|2 



(39) 



(40) 



(41) 



d n denoting (h ■ V). Here it should be remarked that £ is 
a parameter, which depends on the spatial structure of 
the mode at hand. Likewise, one may compute 1%: 



^ = —2 ffda||V t ^°)|| 2 



(42) 



Now using the 2D Green's theorem together with equa- 
tion J3EJ), one gets: 



h = 2eJJda|V< 
5 



(o)i 2 



and, eventually, 



e / Ld 



tfda H , 



(0)|2 



(43) 



(44) 



Considering that the total electromagnetic energy stored 
in the cavity is given by: 



H/ = |JJJd,||^!| 2 = yJJda|^)| 2 
v 5 

the FWHM of the resonances is given by 



(45) 



(46) 



In our context fi c / (i is practically equal to unity. Thus 
one finally has to consider two distinct types of ohmic 
losses: those located at the surface of both ends, which 
amount to attenuation along propagation, and ohmic 
losses upon reflection at the contour: 



r 



ohm 



pprop , p 



rcfi 



(47) 



pprop = 2 / LP _ S cnds U 

d V 2^,a ends d 



^rcfi 



Lo cont u> 
2A 



(48) 



(49) 



Here, we have introduced two different effective conduc- 
tivities (<7ends and <7 con t) and their corresponding skin 
depths (<5 en( j s and <5 cont ) to account for two different types 
of copper used in our experiment, and for possible differ- 
ent surface states for the top and bottom plates and the 
inner surface of the copper frame used as the contour. 
These two contributions to the widths are also quite dis- 
tinct in their physical interpretation. Indeed, r prop truly 
corresponds to losses endured by a plane wave propagat- 
ing in free space between two parallel infinite metallic 
planes. It is a slowly varying function of frequency, de- 
pending neither on the transverse geometry of the cavity, 
nor on the spatial distribution of the wavefunction in the 
cavity. r rcfl , to the opposite, is related to a loss mecha- 
nism located on the contour of the cavity, which clearly 
depends on the geometry of the latter, and chiefly, on the 
spatial distribution of the normal gradient of the wave- 
function via the quantity £. r rcfl therefore fluctuates from 
mode to mode and, in the case of the rectangular cavity, 
its explicit form will be given in the following section. 
By using a boundary perturbation technique to compute 
losses pertaining to reflections on the contour it is shown 
in appendix iBl that, to each correction of the imaginary 
part of the frequency due to ohmic losses, there is a cor- 
respondingly equal correction of the real part. Moreover 
this boundary approach sheds light on the intimate con- 
nection between the boundary conditions on the contour, 
leading to non-purely real wavefunctions, and the fluctu- 
ating partial widths r refl plT|. 

Collecting the above results with those obtained in the 
previous section we are now in a position to write the 
S'-matrix between weakly coupled pointlike antennas in 
a 2D cavity in the presence of ohmic losses. It reads: 



Sab = 5 ab - 2iT 2 H J2 



i w 2-c4 + i4 0) (r$? m + rr) 

(50) 

where w„ = uj^ — r° hm /2 and the }'s are the unper- 
turbed eigcnfrcqucncics of the ideal lossless cavity, and 
where 



T{uj) = c 



Zq sin 2 ^/ 

V ojZ* „ rtP , u)l* 7 



2c 



(51) 



the contribution of antennas to the widths being given, 
at leading order, by 



T 2 (uj„) 



El* 



(0), 



(0) A^^n VCJ 



(52) 
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Here it should be remarked that the second factor ip in 
equation (|50|) should not be a complex conjugate. Indeed, 
due to ohmic losses, the wavef unctions are no longer real 
and the S-matrix cannot keep its unitarity. Nonetheless, 
it obviously has to remain symmetric. In the previous 
section we used the Hermitian formalism for the sake 
of convenience, but it turns out to be inappropriate for 
the present purpose (see for instance reference 13] about 
self-adjoint systems). 

IV. EXPERIMENTAL VALIDATION IN A 
RECTANGULAR CAVITY 

A. A preliminary global test 

The aim of this section is to check the pertinence of 
the description given above in a rectangular cavity where 
cigenfunctions and cigcnfrcqucncies arc easily calculated 
in the limit of vanishing losses. A preliminary test, which 
does not involve any sophisticated fitting procedure, con- 
sists in comparing the average transmission between two 
antennas with the corresponding quantity deduced from 
equation (|50|l . According to this equation, the transmis- 
sion coefficient (a 7^ b) for lu = u n approximately reads: 

I c /, , \| ^ n T 2/, , \ IV'n(^a)^'n(^i>)l frn\ 
\b ab {U! n )\ ~ 21 (W n ) Tqj , (53) 

where T n = r° hm + r^ nt . In a lossless rectangular cavity 
with sides L x x L y , the eigenfrequencies are: 

where Z, m are integers. The corresponding eigenfunc- 
tions read: 

^i%( x i y) = -j== sin sin . (55) 
Assuming that |-0„(r)| ~ ^^(x, y)\, one deduces: 

T*(u n )^ LxL ^ nTn (\S ab (u; n )\), (56) 

where the average is performed on the positions f a and 
f% of the antennas. In our experiments the average was 
obtained by performing the 45 distinct transmission mea- 
surements that the 10 antennas allow. The values of w„ 
and r„ used in the experimental evaluation of the right- 
hand side of 1|56[) were obtained, in a rough way, through 
the analysis of the group delay. This quantity, defined as 
the derivative of the phase ip a b of S a b with respect to uj, 
presents rather well defined extrema at frequencies close 
to eigenfrequencies. For a gross estimation of T n we used 
r„ ~ 2/(diy9 a h(a; n )/dw), which is exact for an isolated 



resonance. Note that, in the case of moderate or large 
modal overlap, this method generally overestimates the 
widths. Knowing the dimensions L x and L y of the cav- 
ity, the right-hand side of l|56|) only depends on the length 
parameters Z e // and I*. Recall that I* is the length of the 
perturbed region at the end of the coax, and that l e f / is 
the effective length over which the field inside the cavity 
is coupled to the antenna. One may assume, in this pre- 
liminary global test, that Z e // remains close to the length 
I of the part of the antenna which lies inside the cavity: 
l e ft~l = 2.0 ± 0.3 mm. In the same way, a rough esti- 
mate of I* is given by the distance between the end of the 
antenna and the reference (calibration) plane of the coax: 
16.5±0.2mm. Figure|3compares T 2 as obtained through 
equation (|51|l with its experimental value deduced from 
equation (|56[l for the first 348 resonances up to 5.5 GHz. 
The lowest and the highest continuous curves correspond 
to l e ff = 1.7mm and l e ff = 2.3mm respectively. At the 
resolution of the presented figure, the curves obtained for 
values of I* ranging from 16.3 mm to 16.7 mm are not dis- 
tinguishable. A fair agreement is observed between the 
experiment and our model. A first consequence of this 
test is to substantiate the correspondence between the 
length parameters of our model and the actual lengths 
of the coax antennas. Note also that this global prelim- 
inary test requires no sophisticated data processing of 
the individual resonances. In the following, we will see 
how our model for coupling with antennas remains quite 
satisfactory when put to more stringent tests. 




Frequency in GHz 



FIG. 3: Experimentally based estimation of T 2 (uj„) as given 
by 156|l (crosses). The gray region corresponds to its theo- 
retical expression l|51|l for values of l e ff between 1.7 mm and 
2.3mm (1* = 1.65cm). 



B. A test with individual resonances 

In the global test presented above we had no need of 
a precise knowledge of the wavefunctions at the anten- 
nas. To check the validity of formula 150JI. we have de- 
veloped an original fitting procedure (see |22|) to extract 
the actual complex eigenfrequencies tu n and the complex 
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amplitudes 

A ab = -2iT 2 (Lu n )iP n {r a )^ n (r b ) . (57) 

With this procedure, based upon a mixture of robust al- 
gorithms, we could check that the actual eigenfrequencies 
remain very close to the unperturbed values given by H54|) 
in the frequency range studied here. Now we proceed to 
verify that the amplitudes A n deduced from our measure- 
ments are well described by ()57|l with expression Q51|l for 
T and formula l|55|l to approximate the true eigenfunc- 
tions. Indeed, even if the existence of non-uniform losses 
(chiefly those associated to r rofl ) leads to non-purely real 
wavefunctions, we will see below that the corrections re- 
main very small. From Q57[l it is easily deduced for three 
different antennas a, b and c: 

\Aab\\ Aacl 

]a\\ = 2T2 M\iPn(r a )\ 2 . (58) 

As there are 36 different ways of combining the 45 am- 
plitudes A ab , yielding 36 slightly different values of (I58|l . 
one can use the following average estimate: 

2T 2 K)|^(x a ,, a )| 2 = JL £ , (59) 

b.c^a ' n ' 
b<c 

where |^>„(f)| ~ iWm^y)! m assumed. For a given 
resonance one can directly compare the above quan- 
tity with 2T 2 (uj n )\%l) < f^ n (x a , y a )\ 2 for the ten antennas 
(a = 1, ...,10). This comparison is shown on figure 01 
for two distinct resonances, namely n = 78 at 2.655 GHz 
and n = 238 at 4.558 GHz. As the great majority of res- 
onances that we are concerned with are narrow enough 
to ensure a good correspondence with the unperturbed 
wavefunctions, we observe a fairly good agreement. 

To extend our comparison to all resonances, one can 
now average over all ten locations of antenna a. For 
the best values of l e ff = 1.9 mm and I* = 16.5 mm ob- 
tained in the frequency range from 2 GHz to 5.5 GHz, 
figure El shows the comparison between the experimental 
quantity (± £ h<c#a \A a n b \ \A?\/\A%\) a and the predic- 
tion 2T 2 {u n )(\^°l(x a ,y a )\ 2 ) a up to 5.5GHz. The agree- 
ment is excellent on the average. By a close inspection, 
e.g. between 2 GHz and 3 GHz as shown in the inset, one 
can notice that the agreement is generally excellent even 
at the level of individual resonances. Rare important dis- 
crepancies are observed for very close neighboring eigen- 
frequencies (quasi-degeneracies) due to modal overlap- 
ping. Indeed, when the latter effect is not negligible, one 
expects that the spatial distribution of the wavefunction 
results from a linear combination of neighboring unper- 
turbed wavefunctions [23|. Beyond 3 GHz this effect de- 
teriorates the agreement due to the concomitant increase 
of the total width and decrease of the mean spacing lead- 
ing to an inadequacy of the zeroth order eigenfunctions 
we use for our test. 



n = 78 


4- 


2T 2 K)|V(r a )| 2 

/A 






- n = 238 




+\ 



123456789 10 



antenna's index 



FIG. 4: Comparison of £23 with 2T 2 (tj n )|^ i (0 j 1 (a;a, y a )\ 2 at 
the ten antennas (indices are introduced in figure for reso- 
nances n = 78 and n = 238. 




Frequency in GHz 



FIG. 5: Comparison between the experimental quantity 
(ssEkc*,!^ hm/\A b n c \) a (crosses) with the prediction 
2T 2 (Lo n )(\ipl°^(x a ,y a )\ 2 } a (line). Inset : enlarged view within 
the range 2-3 GHz. 

C. Partial width decomposition 

As seen above, the total width of a resonance can be 
decomposed as a sum of three partial widths associated to 
losses through the antennas, ohmic losses on the contour 
of the cavity, and ohmic losses at the surface of both 
ends which appear as damping along propagation. In a 
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rectangular cavity, the total width of the nth resonance, 
characterized by the quantum numbers I and to, is given 
by 

r„ = r^ + r;-«+rp ro PK). (60) 

By using expression i|55|) for the wavefunctions in order 
to evaluate the factor £ in equation (|41|) . one obtains 

Ac 2 ft 2 m 2 \ , . 

whence, using 

™a c 2 <5 ron t(w„) / I 2 to 2 \ , . 

These partial widths clearly vary from mode to mode as 
illustrated in figure where — 1 is shown for each 
cigcnfrcquency up to 5.5 GHz. Note that the £z, m 's oscil- 
late around unity and vary at most by 23 %. 



u3 




1 2 3 4 5 6 
Frequency in GHz 

FIG. 6: Values of (£;. m — 1) after equation 16111 for all reso- 
nances up to 5.5 GHz. 

The widths associated to losses through antennas also 
vary from mode to mode: 

r^ = ^£l<tol 2 - (63) 

W i,m c=l 

In the previous subsection we already checked the va- 
lidity of the above formula for all resonances shown on 
figure 13 since Tf 1 ^ is essentially proportional to the sum 
of expression l|59fl over all antennas. 

Thus by fitting the experimental transmission by for- 
mula H5Ufl . one obtains a direct measure of the total width 
and an indirect measure of T™* , thus enabling us to eval- 
uate the two effective conductivities (Tends and <7 CO nt. A 
representation of all the ohmic widths r° hm = r„ — Tf"' 
up to 3 GHz is given in figure [7| A comparison is shown 
between theoretical ohmic widths (crosses) and experi- 
mental ohmic widths (continuous curve). The smooth 
dominant contribution given by r prop (a;) is indicated by 
the dashed curve. The agreement on the mean level and 



on the amplitude of fluctuations (only present in r™ fl ) is 
fairly good whereas, resonance by resonance, the agree- 
ment is not systematic, essentially due to the effect of 
modal overlap mentioned above. 



1 




'0 0.5 1 1.5 2 2.5 3 

Frequency in GHz 

FIG. 7: Behavior of the ohmic widths up to 3 GHz. Compar- 
ison is shown between theoretical ohmic widths (crosses) and 
experimental ohmic widths (continuous curve). The dashed 
curve indicates the smooth contribution r prop (a>). 



V. CONCLUSION 

In conclusion, we have tried to provide a better under- 
standing of the physical mechanisms at the origin of reso- 
nance broadening in microwave cavities. We have explic- 
itly developed an S-matrix model including the frequency 
dependent coupling of the antennas and accounting for 
ohmic absorption at the boundaries of a two-dimensional 
cavity. We have especially emphasized the necessity to 
distinguish between ohmic attenuation along propaga- 
tion, leading to a smooth frequency dependent contri- 
bution to the total width, and localized absorption at 
the contour of the cavity, yielding a contribution varying 
from mode to mode. We have performed experiments 
where we analyzed the transmission versus frequency in 
terms of a Brcit-Wigncr decomposition deduced from our 
model. In the rectangular cavity we used, all the quanti- 
ties involved in our theoretical description could easily be 
calculated in the perturbative limit of small or moderate 
modal overlap. We therefore have been able to validate 
our approach and provide a very precise estimation of the 
various contribution to the total widths. In particular, 
the varying contribution of ohmic losses at the contour 
could be quantitatively checked at the level of individual 
resonances, except for quasi-degenerate modes. This ap- 
proach has recently enabled us to relate the losses at the 
contour, in a chaotic cavity, to the imaginary part of the 
wavef unct ion [2lj . 

We believe that our present approach is an important 
step to test existing or yet to come theories of open or ab- 
sorptive chaotic wave systems. Indeed, such theories gen- 
erally assume that losses are associated to distinct well 
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identified coupling channels |ll|. It is nonetheless not ob- 
vious that these can be used to describe different sources 
of loss as damping along propagation or ohmic dissipa- 
tion at the contour. For instance, absorbing boundaries 
may be viewed BjS cl number (of the order of L/X) of 
distributed coupling channels. Indeed, to mimic absorp- 
tion, recent theoretical predictions have been proposed 
but only in the asymptotic limit of a lar ge n umber of ef- 
fective channels with vanishing couplingpil |25| (see also 
0). 
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The eigenfunctions f a (f) obey boundary Dirichlet con- 
ditions, and the last term in (|A4|) vanishes for z — » oo. 
The self-adjointness condition can thus be written: 

M 

^[(^(0)^(0)-^(0)^,(0))+... 

c=l 

+ ip^iOm (k)ip al (f c ) - T c (fc)<(r c ) m (0)] =0. 

(A5) 

This condition is non-trivially fulfilled by imposing: 



with P. Sebbah, R. L. Weaver and H.-J. Stockmann. 

APPENDIX A: EVALUATION OF THE 
5-MATRIX 

1. Self-adjointness condition 

In the Hilbert space of the complete problem (cavity 
and antennas), one defines the following scalar product: 

(Al) 
(A2) 

(A3) 
(A4) 



2. The coupling matrix elements 

Assuming Eq = fc 2 , equation (|20|) now reads: 

M 

(A + fc 2 V Q (r) =^T c (fe)J(f- ? )(p P M ■ (A7) 

c=l 

Furthermore, the Green's functions of the isolated cavity 
are given by: 



With this scalar product, the self-adjointness condition, 
reads 

(«*,#') - (§,«*') = JJ df [^*(r)(A^(r)) - (Aip a (r)ycp a ,(r)] + 

M M 



+ 53(T c *(*)^ o (0)^(f c ) - <Pl iro)T c {k) m {Q)) -J2 dz c [<pZ(z c ) m (z c ) - ^> c )^,(z c )] 



c=l c=l ' 

Using Green's theorem for the first term and integrating the third one by parts, one obtains: 

(W*,* 7 ) - (*,«#') = jd.ft[(p a {f)iy<pa,{f)) - (V(^ a (r))^(r)] +... 
M M 



(A6) 



(A + k 2 )G(r 7 r',k 2 ) = 5{r - r') . (A8) 
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By comparing equations 1|A7() and (|A8|) . for f — f b , one 
obtains: 



M 



,{r b ) =J2 T c(k)G(n,f c ,k)^ c (0) , (A9) 



With the boundary conditions (|A6|I and the expression 
(|15fl of tpp, relation (|A9() becomes: 



m . M 

A c m - l T T* b (k) J2 T c (k)G(f b) f c , k)Al = A c out + jT b *(k) J2 T c (k)G(n, r c , k)A 



out J 



r 



6=1 



(A10) 



or, in a matrix form. 

An - iKAin = Aout + iKAout , (All) 

where K is an M x M matrix defined by: 



K ab (k) = ^G(r a ,r b ,k) Tb{k) 



Vk 



Vk 



(A12) 



From Aout = SAin one thus makes the relation l|26(l 
between K and the S'-matrix explicit. By introducing in 
(|A12|) the expansion of the Green's function in terms of 
the eigenfunctions ip„ (r) = (r\v) of the isolated cavity: 



N 



G(r,r',E)=J2 



k 2 - kl 



(A13) 



one finds the following expression for the i-Gmatrix 

Kab{k) = {-^k—) ^ki {-^k— 

(A14) 

Finally, the elements of the coupling matrix W, related 
to K by equation (|27|l . are given by: 



IF, 



T c (k)rJr c ) 



fiC. 



(A15) 



APPENDIX B: PERTURB ATIVE BOUNDARY 
CONDITIONS 

An alternative way of computing the losses at the con- 
tour is easily obtained by following a boundary perturba- 
tion technique (see for instance 20]). Indeed, calling ip 



the solution at the real eigenfrcquency uio of the zeroth- 
order problem defined by equations and l|37|) . the 
perturbation of the boundary conditions can be written: 



ip ~ -(1 + i)^d n i> on C . 



(Bl) 



Green's theorem applied to this 2D problem straightfor- 
wardly yields: 



-(i + O 



1 fi c S § c d£ \d n tp \ 



(B2) 



eii 2ii SS s da\tP \ 2 ' 
Writing the perturbed eigenfrequency as cj = Wo + Su> — 
iT/2, equation (|B2|) becomes: 



Scj-iT/2~ -(1 + i) 



1 {i c 5 § c M\d n %\ 
efiA^iuJo JJ 5 da|i/'o| 2 



(B3) 



leading, in the case at hand, to equal corrections on both 
real and imaginary parts of u>. This last comment is gen- 
eral and applies as well to the corrections originating from 
ohmic losses at top and bottom of the cavity. Thus, this 
pcrturbative approach completely agrees with the one de- 
veloped in section lTTll as long as the widths are concerned 
and completes it as it provides an estimation of the fre- 
quency shift of the resonances related to ohmic losses. It 
is also particularly interesting to note that equation JB1|) 
gives a quite natural hint of how ohmic losses on the 
contour induce a small amount of complexity for wave- 
functions that are purely real in the unperturbed limit. 
As long as a perturbative approach is valid, an immedi- 
ate connection is deduced between the complex character 
of wavefunctions and the fluctuating part of the widths 
embodied in the quantity £ defined in (|4T1) 2~Q] . 
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